Entanglement oscillation and survival induced by non-Markovian decoherence 

dynamics of entangled squeezed-state 
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We study the exact decoherence dynamics of the entangled squeezed state of two single-mode 
optical fields interacting with two independent and uncorrelated environments. We analyze in detail 
the non-Markovian effects on the entanglement evolution of the initially entangled squeezed state for 
different environmental correlation time scales. We find that the environments have dual actions on 
the system: backaction and dissipation, fn particular, when the environmental correlation time scale 
is comparable to the time scale for significant change in the system, the backaction would counteract 
the dissipative effect, fnterestingly, this results in the survival of some residual entanglement in the 
final steady state. 
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I. INTRODUCTION 

Studies on the decoherence dynamics of open quantum 
systems are of great importance to the field of quantum 
information science Any realistic analysis of quan- 
tum information protocols should take into account the 
decoherence effect of the environment. In many quantum 
communication and computation schemes, information is 
transmitted using photons. For instance, the first ex- 
perimental verification of quantum teleportation Q used 
pairs of polarization entangled photons to transfer the 
polarization state of one photon onto another Q . Within 
a year, unconditional quantum teleportation of optical 
fields was demonstrated experimentally using squeezed- 
state entanglement 4, 5] . Given their central role in these 
schemes and many others, much work has been carried 
out on the decoherence dynamics of optical fields. In 
particular, several authors have studied the continuous 
variable entanglement o f op tical fields (see, for instance, 

Refs. @H Ella). 

Conventional approaches not only treat the interac- 
tions between the quantum system S of interest and its 
environment E perturbatively, they also assume that the 
environmental correlation time te is small compared to 
the time scale To for significant change in S. These yield 
approximate equations of motion, i.e. master equations, 
under the Born-Markov approximation [l3|, [3] • Indeed, 
many studies on the entanglement dynamics of continu- 
ous variable system relied on this approximation @, 0, Q . 
However, it is evident from recent experiments (see, for 
instance, Refs. |15l Il6l. Il7|). that there are many phys- 
ically relevant situations where the Markovian assump- 
tion does not hold, and a non-Markovian treatment of 
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the open system dynamics is necessary. So, there has 
been an increasing interest in the understanding of the 
decoherence effect of open quantum system going beyond 
the Born-Markovian approximation in the last decades 

Very recently, some phenomenological models on non- 
Markovian entanglement dynamics of optical fields have 
been investigated |9l. IToL ITTj| . It was found that in contrast 
to the monotonic decrease of entanglement over time in 
Born-Markovian entanglement dynamics @, 0, there 
are transient entanglement oscillations in non-Markovian 
ones. These oscillations are caused by the backactions 
of the environments on their respective local quantum 
systems [l^, [H| • The backaction, characteristic of non- 
Markovian dynamics, means that the environments with 
their states changed due to interactions with the systems, 
in turn, exert their dynamical influences back on the sys- 
tems. 

In this paper we consider the exact decoherence dy- 
namics of the continuous variable entangled squeezed 
state of two single- mode optical fields, Si and S2, that 
are spatially separated. Each optical field, Sk, interacts 
with its own environment Ek (k = 1, 2). E± and E2 are 
independent and uncorrelated. We study the exact en- 
tanglement dynamics of the two optical-field system for 
different te's in comparsion with to, and analyze when 
the system dynamics will exhibit novel non-Markovian ef- 
fects, and provide a detailed description of these. To this 
end, we use the influence functional formalism [l9|, [20| . 
developed explicitly in Refs. [H, [22j. Our results 
show that besides the short-time oscillations, the non- 
Markovian effect can affect the long-time behavior of the 
system dynamics and the steady state as well. In partic- 
ular, when te is comparable to To, we find that the back- 
action effects counteract the dissipative effects of E\Ei 
on S1S2 respectively. This leads to there being some 
nonzero residual entanglement in the steady state. 

Our paper is organized as follows. In Sec. II, we intro- 
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duce a model of two single-mode optical fields in two 
independent and uncorrelated environments; and out- 
line the exact dynamics that was derived in detail in 
Ref. [22]. In Sec. Ill, using logarithmic negativity as 
an entanglement measure of continuous- variable states, 
we discuss the entanglement dynamics of the entangled 
squeezed state. Sec. IV presents the numerical results of 
the entanglement dynamics, where we analyze explicitly 
the non-Markovian effect of the environments on the sys- 
tem for different r^'s in comparison with t . Finally, we 
conclude in Sec. V. 



II. THE TOTAL HAMILTONIAN AND EXACT 
REDUCED SYSTEM DECOHERENCE 
DYNAMICS 

The total Hamiltonian of the system S1S2 plus envi- 
ronment E1E2 is given by 

H = H S + H E + H l , (1) 

where 

2 

H s = fiw k a\a k , 

k=l 
2 

H E = J2^2 tiWklb l:l bkl ' 

k=l I 
2 

Hi = J2J2 h ^kiatb kl +g* kl a k bl), (2) 
fc=i 1 

are, respectively, the Hamiltonian of the two optical 
fields, the two independent environments, and the in- 
teractions between them. The operators a k and a\ (k = 
1, 2) are respectively the annihilation and creation oper- 
ators of the fc-th optical mode with frequency u> k . The 
two independent environments are modeled, as usual, by 
two sets of harmonic oscillators described by the annihi- 
lation and creation operators b k i and b kl . The coupling 
constants between the fc-th optical field and its environ- 
ment are given by g k i. Currently, most quantum optical 
experiments are performed at low temperatures and un- 
der vacuum condition. In this case, vacuum fluctuations 
are the main source of decoherence. Therefore, we take 
the environments to be at zero temperature throughout 
this paper. 

Since we are only interested in the dynamics of S1S2, 
we like to eliminate the degrees of freedom of E\E2- The 
influence-functional theory of Feynman and Vernon (l9j 
enables us to do that exactly. By expressing the for- 
ward and backward evolution operators of the density 
matrix of the system S1S2 plus environment E1E2 as 
a double path integral in the coherent-state representa- 
tion [23l | , and performing the integration over the degrees 
of freedom of E1E2, we incorporate all the environmen- 
tal effects on S1S2 in a functional integral named influ- 
ence functional [H, 0, [HJ . The reduced density matrix, 



which fully describes the dynamics of S1S2 is given by 

p(a.f,a!f)t) = J dp{a i )dp{a' i )J(af,a.' f ]t\a. t ,a.' t ]Q) 
X /j(ai,c^;0), (3) 

where p(a.f,<y!j]t) = {af\p{t)\a!^) is the reduced density 
matrix expressed in coherent-state representation and 
J(a.f, a!f] t\ot.i, a^; 0) is the propagating function. In the 
derivation of Eq. (O, we have used the coherent-state 
representation 

2 

|a) = \\ \a k ), \a k ) = exp(a fc aj £ )|0 fe ). (4) 
fc=i 

which are the eigenstates of annihilation operators, i.e. 
<ik\ctk) = ct k \a k ) and obey the resolution of identity, 
J dp (a) |q;) (a I = 1 with the integration measures de- 
fined as dp. (a) = J]; e ~ a '°" ds ^ 1 . a denotes the com- 
plex conjugate of a. 

The time evolution of the reduced density ma- 
trix is determined by the propagating function 
J7(a/,aj;t|aj,a^;0). The propagating function is ex- 
pressed as the path integral governed by an effective ac- 
tion which consists of the free actions of the forward and 
backward propagators of the optical-field system and the 
influence functional obtained from the integration of en- 
vironmental degrees of freedom. After evaluation of the 
path integral, the final form of the propagating function 
is obtained as follows 

2 

J(af,a' f ;t\a i ,a' i ;0) = exp j ^ [u k (t)a kf a ki 

k=l 

+ u k {t)a' kt a' kf + [1 — \u k (t)\ 2 }a ki a ki ] j, (5) 
where u k {r) satisfies 

u k (r) + iuj k u k (r) +/ p k (r - r')u k (r') = (6) 
Jo 

with p k {x) = J2i e~ lu ' x \g k i\ 2 being used. Combining Eq. 
([5]), we can get the exact time-dependent state from any 
initial state by the evaluation of the integration in Eq. 

©■ 

To compare with the conventional master equation de- 
scription of such system, we now derive a master equation 
from the above results. After taking the time derivative 
to Eq. ([3]) and recalling the explicit form of Eq. (O, we 
can derive an exact master equation 

2 

p(t) = ~[H'(t),p(t)]+J2^k(t)[^kP(t)al 

-a\a k p(t) - p{t)a\a k ], (7) 

where 

2 

H'(t)=Y,1& k (t)a\a k , (8) 
fc=i 
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is the modified Hamiltonian of the two optical modes and 



u k (t) 



-T k {t)-iVt k {t). 



(9) 



Eq. ([7]) is the exact master equation for the optical-held 
system, fl k (t) plays the role of a time-dependent shifted 
frequency of the fc-th optical held. Tf, (t) represents the 
corresponding time-dependent decay rate of the held. We 
emphasize that the derivation of the master equation goes 
beyond the Born-Markovian approximation and contains 
all the backactions between the system and the environ- 
ments self-consistently. All the non-Markovian character 
resides in the time-dependent coefficients of the exact 
master equation. 

The time-dependent coefficients in the exact master 
equation, determined by Eq. (J9j) , essentially depend on 
the so-called spectral density, which characterizes the 
coupling strength of the environment to the system with 
respect to the frequencies of the environment. It is de- 
fined as Ji{lo) — J2 k \gik\ 2 5(u! — coi). In the continuum 
limit the spectral density may have the form 



Jfe(w) = rj k cu( — ) 



ra-1 



(10) 



where lo c is an exponential cutoff frequency, and r\ k is 
a dimensionless coupling constant between S k and E k . 
The environment is classified as Ohmic if n — 1, sub- 
Ohmic if < n < 1, and super-Ohmic for n > 1 [2(| 
|24 |. Different spectral densities manifest different non- 
Markovian decoherence dynamics. 

We note that our exact master equation reduces to the 
conventional master equation under the relevant Markov 
approximation. The coefficients in the master equation 
([t|) become time-independent [HJ 



T k (t) = TrJ k (u>k), 

n*(t) = Lo k -v i 



J(uj)duj 



J Q uj — u>k 

where V denotes the Cauchy principal value. The coeffi- 
cients in Eqs. (jTTJ) are precisely the corresponding ones 
in the Markovian master equation of the optical system 

m. 



III. THE NON-MARKOVIAN 
ENTANGLEMENT DYNAMICS OF THE 
ENTANGLED SQUEEZED STATE 

Initially at time t — 0, S1S2 is in an entangled squeezed 
state. The entangled two-mode squeezed state is defined 
as the vacuum state acted on by the two-mode squeezing 
operator 

|V>(G)) = exp[r( ai a 2 - <4<4)]|00), (12) 

where r is the squeezing parameter. In the coherent-state 
representation, this initial state is given by 

exp[— tanhr(aiia 2 i + Q.' u a' 2i )] 



p(a i; a-;0) 



cosh 2 r 



(13) 



The state approaches the ideal Einstein-Podolsky-Rosen 
(EPR) state [25[ in the limit of infinite squeezing (r — > 
00). The traditional way to generate the entangled two- 
mode squeezed state is via the nonlinear optical pro- 
cess of parametric down-conversion [261] . Recently, a 
microwave cavity QED-based scheme to generate such 
states has also been proposed (27|. After generating the 
entangled state given by Eq. (fT2"|) , the two cavity helds 
are then propagated, respectively, to the two locations 
separated between the sender and the receiver. A quan- 
tum channel is thus established through the entangled 
two-mode squeezed state and is ready for teleporting un- 
known optical coherent states 0, H[ . 

At t > 0, due to interactions with E1E2, \i>(0)) evolves 
to a mixed state. A straightforward way to obtain the 
time-dependent mixed state is by integrating the prop- 
agating function over the initial state of Eq. ([3]). Then 
the time-evolution solution of the reduced density matrix 
can be obtained exactly, 

p(a f ,a' f ;t) = aexp[^ (^a k fa k 'f+ca k fa' kf +^-a' kf a' k/f )], 



k=tk' 



where 



cosh 2 |r| [1 - tanh 2 \r\ (1 - \u(t)\y] 



— tanh |r| u(t) 2 



1 -tanh 2 \r\ (1 - \u{t)\ z f 
tanh 2 \r\ (1 - |w(i)| 2 ) \u{t)\ 2 
1 - tanh 2 |r| (1 - Ki)| 2 ) 2 ' 



(14) 

, (15) 
(16) 

(17) 



To measure the entanglement in continuous variable 
system, one generally uses the logarithmic negativity [28[ . 
The logarithmic negativity of a bipartite system was in- 
troduced originally as 



E 



N 



log. 



2 E 



a; 



(18) 



where A~ is the negative eigenvalue of p Ti , and p Ti is a 
partial transpose of the bipartite state p with respect to 
the degrees of freedom of the i-th party. This measure 
is based on the Peres-Horodecki criterion [29|, [3(| that 
a bipartite quantum state is separable if and only if its 
partially transposed state is still positive. 

For the continuous-variable (Gaussian-type) bipartite 
state, its density matrix is characterized by the co- 
variance matrix defined as the second moments of the 
quadrature vector X — (xi,pi, ^2,^2), 



Vu 



(AXiAXj + AXjAXi) 



(19) 



where 



AXi = X, - (Xi), and x { = P 



1V2 



The canonical commutation relations take the form as 
[X h Xj] = iUij, with U = I Q j\ and J = I j J 
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defining the symplectic structure of the system. The 
property of the covariance matrix V is fully determined 
by its symplectic spectrum v = {y\, V2), with ±i/j > 0) 
the eigenvalues of the matrix: iUV. The uncertainty 
principle exerts a constraint on Vi such that i/j i [3 11 ] - 
Thus the Peres-Horodecki criterion for the continuous- 
variable state can be rephrased as the state being sepa- 
rable if and only if the uncertainty principle, V + ^U ^ 0, 
is still obeyed by the covariance matrix under the par- 
tial transposition with respect to the degrees of freedom 
of a specific subsystem [32} . In terms of phase space, 
the action of partial transposition amounts to a mirror 
reflection with respect to one of the canonical variables 
of the related subsystem. For instance, V = AVA, and 
A = diag(l, 1, 1, —1) is the partial transposition with re- 
spect to the second subsystem. If a Gaussian-type bipar- 
tite state is nonseparable, the covariance matrix V will 
violate the uncertainty principle and its symplectic spec- 
trum v — V2) will fail to satisfy the constraint Vi ^ ^. 
The logarithmic negativity is then used to quantify this 
violation as 12811 





FIG. 1: Damping rate T(t) and logarithmic negativity -Bjv(t) 
for the entangled squeezed state as function of dimension- 
less quantity tu^t and their corresponding Markovian results 
(dashed lines). The parameters u) c /uq = 50.0, r\ = 0.1, and 
r = 1.0 are used in the numerical calculation. 





E N = max{0, - log 2 (2j> min )}, 



(20) 



where j> m i n is the smaller one of the two symplectic eigen- 
values It is evident from Eq. ([20]) that, if V obeys 
the uncertainty principle, i.e., ^ i, then E^{p) = 0, 
namely, the state is separable. Otherwise, it is entan- 
gled. Therefore, the symplectic eigenvalue //,„;„ encodes 
a qualitative feature of the entanglement for an arbitrary 
continuous- variable bipartite state. 

With this entanglement measure at hand, we study 
now the entanglement dynamics of the squeezed-state 
quantum channel in our model. From the time-dependent 
state, the covariance matrix for the optical field can be 
calculated straightforwardly, 



V 



a(i+Q 

2(l-d) 2 



qR.c[fr] 





1/(1+0. 

2(l-d) 2 
alm[6] 



aim [b] — aRc[6] 



aRc[b] 

x 

aIm[Z>] 

i/(i+0 

2(l-d) 2 





aim [b] 

x 

-aRc[b] 

x 



2/(1+0 
2(l-d) 2 
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, and d = c 



-^1 



(21) 



And 



where x — [(1— c) 2 — \b 
the logarithmic negativity Ejsr(t) can also be obtained 
exactly from Eq. ([2TJ]) . It is easy to verify that the initial 
entanglement is £W(0) = 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

In the following, we analyze explicitly the exact deco- 
herence dynamics of the entangled squeezed state of S1S2 
under the influence of E\E2- For simplicity, we assume 
from here on that the two optical fields are identical, i.e., 
u>i = L02 = ^o! an< i they interact with the same strength, 
.91/ = 52/ = 9u with their individual environments. For 



FIG. 2: Damping rate T(t) and logarithmic negativity -Bjv(t) 
for the entangled squeezed state as function of dimension- 
less quantity uot and their corresponding Markovian results 
(dashed lines). The parameters lj c /ujq = 1.0, r\ = 5.0, and 
r = 1.0 are used in the numerical calculation. 



definiteness, we consider both E\ and E2 to have Ohmic 
spectral density. The environmental correlation time te 
in this case is roughly inversely proportional to the cut- 
off frequency cu c in Eq. (flO|) . i.e., te — l/w c [18]. It is 
emphasized that the cutoff frequency u> c , which is orig- 
inally introduced to eliminate infinities in frequency in- 
tegrations, therefore also determines if the dynamics of 
open system 5* is Markovian or non-Markovian. Our non- 
perturbatively derived exact results allow us to explore 
all these possibilities. 

In Fig. [TJ we plot the numerical results of the de- 
cay rate T(t) and logarithmic negativity £jv(t) when 
te <C to. The positivity of T(t) throughout the whole 
evolution process guarantees the monotonic decrease of 
Sjv(t). Accordingly, the entangled squeezed state even- 




FIG. 3: Damping rate T(t) and logarithmic negativity -Bjv(t) 
for the entangled squeezed state as function of dimension- 
less quantity tu^t and their corresponding Markovian results 
(dashed lines). The parameters w c /wo = 0.2, r\ = 5.0, and 
r = 1.0 are used in the numerical calculation. 
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tually evolves to a product state, namely the ground 
state of the system: p g = |00)(00|. Clearly, in this case, 
the backactions of E1E2 have a negligible effect on the 
dynamics of S1S2, and we say the system dynamics is 
mainly governed by the dissipative effect of the environ- 
ments. There is thus no qualitative difference between 
the exact entanglement dynamics and the Markovian re- 
sults. Quantitatively, however, we note that for t < te, 
the distinctive increase of T(t) results in E^if) decreas- 
ing rapidly. This non-Markovian effect only shows up in 
a very short time scale. In fact, for t > te, T(t) de- 
creases and approaches gradually to a constant value as 
t approaches tq; and the rate of decrease of .Ejv(t) de- 
creases. 

Fig. [2] shows r(t) and E^(t) when te = To- In this 
case, the backactions of E\Ei have a considerable im- 
pact on the dynamics of S1S2, and the Markovian ap- 
proximation is not applicable. Firstly, we note that T(t) 
can take negative values. Physically, this corresponds to 
the systems reabsorbing photons from the environments, 
which leads to an increase in the photon number of the 
systems [lij]. These negative decay rates provide evi- 
dences for backactions in non-Markovian dynamics [34| . 
Secondly, we observe that T(t) approaches zero asymp- 
totically. Both results clearly differ from the Markovian 
ones. Consequently, E^{t) presnts distinctive behaviors 
that are absent in the Markovian results. Firstly, due to 
the negative decay rates, E^{t) shows oscillations. We 
must emphasize that these oscillations are fundamentally 
different from the transient entanglement oscillations pre- 
viously obtained when two optical-fields interact with a 
common environment [l2| . They are caused by the back- 
actions of the environments on their respective local opti- 
cal fields and are characteristic of non-Markovian dynam- 
ics. Similar oscillations have been obtained in a system 
of two two-level atoms in two separated damping cavi- 
ties [Hj]. Secondly, and more interestingly, we find that 
there is some residual entanglement left in the steady 
state. From previous studies [3, [Hl[l3]> one would have 
concluded that non-Markovian effects only show up in 
short-time dynamics. Our results, however, clearly show 
on the contrary that non-Markovian effects can also have 
an influence on the long-time behavior of the system dy- 
namics and the final steady state of the system. This 
counterintuitive behavior can be explained with the fact 
that the dissipative influence on the entanglement dy- 
namics by the environments is strongly counteracted by 
the effect due to their backactions. Consequently, the de- 
cay of the entanglement ceases when the system evolves 
to some steady state, which is not the ground state p g . 

The results when te 3> To are shown in Fig. [3j a situa- 
tion considered in Ref. [TJJ . Due to extremely long mem- 
ory of the environments, the backactions on the systems 
are so strong that they govern the decoherence dynamics. 
As a result, T(t) and hence E^{t) oscillate over a very 
long duration. These oscillations persist even as the state 
approaches the ground state. The 'equilibrium' position 
for the oscillation of T(t) is not at zero, but a small pos- 



itive value. This positivity means the systems dynamics 
will experience a weak dissipation, which is verified by 
the time evolution of Em (t) in Fig. [3J 

In summary, we have studied the exact entanglement 
dynamics of the two optical-field system for different 
T^'s in comparison with tq. Specifically, we have ana- 
lyzed when the system dynamics will exhibit novel non- 
Markovian effects, and provided a detailed description of 
these. 



V. CONCLUSIONS 

We have applied the influence-functional method of 
Feynman and Vernon to investigate the exact entangle- 
ment dynamics of two single- mode optical fields S1S2 
coupled to two independent and uncorrelated environ- 
ments E\Ei- From our analytical and numerical results, 
it is seen that E\Ei exert two competing influences on 
our system. One effect, V, is dissipative and is responsi- 
ble for the decoherence of S1S2. The other, B, is due to 
the backactions of E\Ei on S\Si- The degree of mani- 
festations of T> and B in the dynamics of S1S2 depends 
on te in comparison with To- For te <C to, T> dominates 
and B only gives rise to a transient coherent oscillation of 
S\S2- The state of S1S2 evolves to the ground state p g , 
which is coincident with the Markovian result. If te = tq, 
the near resonant interaction between S1S2 and E1E2 re- 
sults in T> and B being comparable and counteract each 
other. These give rise to transient negative decay rates 
and asymptotical zero decay rate. The state of S1S2 thus 
evolves asymptotically to some steady state, which is not 
the ground state p g . Finally, when te ^> tq, B dominates 
and governs the dynamics of Si £2- The decay rates of 
the system oscillate about some non-negative equilibrium 
position over a very long duration. This slight positiv- 
ity guarantees an overall weak dissipative effect on the 
system dynamics. Therefore, the state of S1S2 eventu- 
ally approaches the ground state with the entanglement 
oscillation persisting on for a very long time. 

The theory we have established is a non-perturbative 
description of the exact decoherence dynamics of a sys- 
tem of two single-mode optical fields. It can serve 
as a useful basic theoretical model in analyzing the 
non-Markovian decoherence dynamics of optical fields 
employed in practical quantum information processing 
schemes. It should be noted that although only the 
Ohmic spectral density is considered here, it is straight- 
forward to generalize our discussion to the non-Ohmic 
cases. 
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